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ABSTRACT 

We present the orbital-phase resolved analysis of an archival FO Aqr observation obtained 
using the X-ray Multi-Mirror Mission (XMM-Newton) , European Photon Imaging Camera (pn 
instrument). We investigate the variation of the spin pulse amplitudes over the orbital period 
in order to account for the effects of orbital motion on spin modulation. The semi-amplitude 
variations are in phase with the orbital modulation, changing from (38. Oil. 8)% at the orbital 
maximum to (13.3±3.7)% at the orbital minimum. The spectral parameters also show changes 
over the orbital period. One of the absorption components increase by a factor of 5 between 
the orbital minimum and maximum. We interpret that this absorption arises from the bulge 
where accretion stream from the secondary impacts the disk. The spectrum extracted from 
the orbital minima and maxima can be fitted with a warm absorber model yielding values Nn 
= 2.09™ x 10 22 and 0.56^°;^ 



x 10 22 cm" 2 ; and log(f) = 0.23±g;l| and <0.30 erg cm s" 1 



respectively, indicating the existence of ionized absorption from the bulge at the impact zone 
which is spread out on the disk. The absorption due to accretion curtain and/or column which 
causes the spin modulation can be distinguished from the disk absorption via spectral modeling. 

Subject headings: binaries:close - Stars: individual: FO Aqr - Intermediate Polars, cataclysmic variables 
- stars:rotation - white dwarfs - X-rays: stars 



1. Introduction 

FO Aqr is a compact binary belonging to the 
sub-category intermediate polars (IPs) of the cat- 
aclysmic variables (CVs) which are composed of 
a white dwarf accreting material from a Roche 
lobe filling main sequence star. IPs have white 
dwarfs with a magnetic field of mostly 1-10 MG 
(although there are a couple of systems with mag- 
netic field strengths up to 30 MG; e.g. Katajainen 
ct al. 2010, Piirola et al. 2008). The accretion 
occurs through a truncated disk and via accretion 
curtains to the magnetic poles of the white dwarf 
(see Patterson 1994; Warner 2003). 

FO Aqr is a well studied IP, with an orbital 
period of 4.85 hr (Osborne & Mukai 1989) and 
white dwarf spin period of 20.9 min (Patterson et 
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al. 1998). It has been observed with almost every 
X-ray mission; EXOSAT, Ginga, ASCA, RXTE, 
XMM-Newton, INTEGRAL, SWIFT and Suzaku 
(Cook, Watson & McHardy 1984, Norton et al. 
1992; Mukai, Ishida k Osborne 1994; Evans et al. 
2004, Parker, Norton & Mukai 2005, Landi et al. 
2009, Yuasa et al. 2010, respectively). Our un- 
derstanding of the accretion scenario in FO Aqr 
has changed over the years: Norton et al. (1992) 
suggested that the system shows diskless accretion 
and that the accretion flow changes poles. How- 
ever, Hellier (1993) argued that there is an accre- 
tion disk in the system and the accretion takes 
place both via a disk and stream overflowing the 
disk. Mukai, Ishida & Osborne (1994) confirmed 
this hybrid accretion mode, but proposed that ac- 
cretion from a partial disk is dominant. Later, it 
was suggested that the accretion mode alternates 
from a hybrid of disk-fed and stream fed accretion 
to disk-fed accretion over the years (i.e. hybrid 
in 1988, disk-fed in 1990 (Norton, Beardmore & 
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Taylor 1996) hybrid in 1993 and 1998 (Beardmore 
et al. 1998), disk-fed in 2001 (Evans et al. 2004)). 
The changes in accretion modes are most likely 
due to changes in the mass accretion rate of the 
system (de Martino et al. 1999). The source shows 
orbital modulations, which are deeper at the lower 
energy regime, indicating that it could be due to 
absorption from structures on the disk and the 
accretion stream (Hellier et al. 1993; Evans et al. 
2004; Parker, Norton & Mukai 2005). The spin 
pulse shape of the system is complicated, with a 
quasi-sinusoiadal component and a notch after a 
"dip" caused by the accretion curtains (Evans et 
al. 2004). The spin pulse profile may vary between 
nearly sinusodial to saw-tooth shapes in different 
observations (Beardmore et al. 1998). 

The X-ray spectrum of the source can be rep- 
resented with a multiple plasma emission compo- 
nent, complex absorption and Gaussian lines (eg. 
Mukai, Ishida & Osborne 1994; Evans et al. 2004; 
Yuasa et al. 2010). A soft blackbody emission 
was not detected by Evans & Hellier (2007) with 
XMM-Newton data or by Yuasa et al. (2010) 
with Suzaku data. However, using the results of a 
combined joint analysis of INTEGRAL/IBIS and 
SWIFT/XRT data, Landi et al. (2009) claimed 
the detection of a blackbody component with a 
temperature of 61 eV which they attributed to the 
irradiation of the white dwarf atmosphere. 

In this work, the spectral and temporal prop- 
erties of FO Aqr over the orbital phase will be 
investigated. In Section 2 observation and data 
preparation will be introduced. In Section 3 the 
spin pulse profile behavior over the orbital phase 
will be outlined. In Section 4 the phase average 
spectrum and the variation of spectral parameters 
over the orbital phase will be investigated. In Sec- 
tion 5 the results of the analyses will be discussed. 

2. Observation and Data Analysis 

FO Aqr was observed for about 35 ks (exposure 
time) on 12 May 2001 (OBS ID: 0009650201) with 
XMM-Newton Observatory (Jansen et al. 2001) 
with all of its instruments on board: a pn CCD 
detector (Striider et al. 2001), two MOS CCD 
detectors (Turner et al. 2001) sensitive in the 0.2- 
15 keV energy range at the focus of three Euro- 
pean Photon Imaging Cameras (EPIC), Reflect- 
ing Grating Spectrometers (RGS; den Herder et 



al. 2001), a high- resolution spectrometer working 
together with the EPIC detectors, and the opti- 
cal monitor (OM) instrument with an optical/UV 
camera (Mason et al. 2001). 

Due to the higher sensitivity compared with the 
MOS, we used the archival EPIC pn data in the 
analysis, obtained using the small window mode, 
with a net count rate of 2.33 ± 0.01 c s _1 in 
the 0.3-10.0 keV interval. Standard data prepa- 
ration and processing was carried out with the 
raw data using the XMM-Newton Science Analy- 
sis Software (sas), version 11.0.0. The extraction 
of the background subtracted source spectrum was 
carried out using especget and the extraction of 
the source and background light curve was carried 
out with evselect which are tools within SAS. 
Data with single- and double-pixel events, i.e. pat- 
terns 0-4 with Flag=0 options were extracted at 
all times. A circular extraction region with a ra- 
dius of 32.4 arcsec was used for the source and 
background photons. To extract the light curve 
and spectra resolved at the orbital phase, we used 
the PHASECALC tool to create a phase column in 
the events file with observation start time being 
the phase. Next, phase-resolved spectra were 
extracted from the modified events file using a 
script. Further analyses of light curve and spec- 
tra were performed with xronos 5.2 and xspec 
12.6.0 packages. We obtained the background 
subtracted source light curve using LCMATH tool in 
FTOOLS v6.10 and applied barycentric correction. 
Our preliminary analysis of FO Aqr is summarized 
in Balman & Pekon (2011). The original analysis 
of the archival XMM-Newton data of FO Aqr can 
be found in Evans et al. (2004). 

3. Spin Modulation Variations over the 
Orbit 

We folded the corrected light curve over the or- 
bital (17458 s) and spin (1254 s) periods (Patter- 
son et al. 1998) using the same epoch given by 
Evans et al. (2004); T =2452041.806 (HJD) for 
the orbital period and T =2452041. 858 (HJD) for 
the spin period. Figure 1 shows the folded light 
curve over the orbital period and the spin period 
separately (see also, Evans et al. 2004). 

In order to investigate the changes in the spin 
pulse over the orbital phase, we extracted light 
curves for each 0.1 orbital phase interval in the 0.3- 



2 




Fig. 1. — The light curve folded over the orbital 
period of 17458 s on top, and the light curve folded 
over the spin period of 1254 s at the bottom. 



10.0 keV range and then folded these over the spin 
period. One orbital-phase bin is 1745.8 s, which is 
greater than the spin period. In general, the pulse 
shapes look similar with spin minima at around 
phase 0.8. However, we find that the amplitude 
of the spin pulses change with orbital phase (See 
Figure 2). To construct Figure 2, we fitted spin 
pulse with a simple sine curve for each orbital 
phase. Then we calculated the percentage vari- 
ations by taking the ratio of the semi-amplitude 
of each fit to the mean value and plotted the per- 
centage variations vs. orbital phase. We also com- 
pared them with the orbital folded light curve. As 
seen in Figure 3, the percentage variation is clearly 
in phase with the orbital motion. The variations 
change from (38. Oil. 8)% at the orbital maximum 
to (13.3±3.7)% at the orbital minimum. 

4. Spectral Variations over the Orbit 

We used the same composite spectral model 
from Evans et al. (2004) in order to fit the phase 
average spectrum of the source. The compos- 
ite model consists of a simple absorption model 
(wabs), two partial covering absorption models 
(pcfabs); three plasma emission models at differ- 
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Fig. 2. — The spin pulse profile of FO Aqr for 
each 0.1 orbital phase interval. The corresponding 
phases are noted on each panel. 




orbital phase 



Fig. 3. — The percentage variation of the spin pulse amplitudes from the mean (below), together with the 
orbital flux variation (above). 



ent temperatures (mekal) and a Gaussian emis- 
sion line at 6.4 kcV (Fe K a ). 

In order to investigate the changes of this com- 
posite model spectrum over the orbital phase, 
we performed orbital phase-resolved spectroscopy. 
We extracted spectra for each 0.1 orbital phase in- 
terval, and fitted them with the composite model 
described above. Table 1 displays the spectral pa- 
rameters over the orbital phases. Figure 4 shows 
the 10 phase resolved spectra extracted from the 
data. The fits to the spectra clearly reveals that 
the neutral Hydrogen column density of the sec- 
ond partial covering absorber PCFABS2 increases 
up to 5 times during the the orbital dipping phases 
(see Figure 5, and Table 1), while the other par- 
tial covering absorber (pcfabsI) shows no partic- 
ular variation (See Table 1). The first two plasma 
temperatures (kTl and kT2) show no significant 
variation over the orbital phase, however the third 
temperature (fcT3) peaks during the orbital min- 
ima indicating spectral hardening (See Figure 5 
and Table 1). 

The shape of the spectrum shows distinct vari- 
ations over the orbital phase as well. During the 
phases where X-ray flux is minimum (i.e. phases 
around 0-0.1) there is a substantial decrease in 



normalized count rates around 1-2 keV regime (see 
Figure 4), while below 1 keV the normalized count 
rates stay at the same level with that of the rest 
of the phases. This soft excess brought our at- 
tention to a possible warm absorber in the line 
of sight. In order to investigate this feature we 
have extracted two spectra from orbital maximum 
and minimum regions (phases between 0.55-0.85 
and 0.9-1.2 respectively), and fitted them with 
the composite model similar to the one in Table 
1, this time replacing the second partial covering 
absorber and one of the plasma emission compo- 
nents with a warm absorber model (i.e. warmabs 
model implemented into xspec). WARMABS mod- 
els the absorption from a photoionized plasma in 
the line of sight with column densities of ions (in- 
cluding small cross sections) coupled through a 
photoionization model using stored level popula- 
tions calculated by XSTAR (Kallman & Bautista 
2001) assuming a given continuum spectrum. The 
fit to the orbital minima and maxima yielded re- 
sults with xt 01 1-39 and 1.24; with values of 
7V H = 2.09±?;^10 22 cm- 2 and 0.561°;*? xlO 22 
cm" 2 ; and log(£) = 0.23t° H and <0.30 erg cm 
s , respectively. Here Nn is the equivalent hy- 
drogen column density of ionized absorption and 
£ = L/n c r 2 is the ionization parameter where 
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Fig. 4. — The orbital-phase resolved spectra of FO Aqr for each 0.1 phase bin interval. The left-hand panel 
shows phases 0.8 (green) 0.9 (dark-blue), 0.0 (red), 0.1 (black) and 0.2 (light-blue). The right-hand panel 
shows the phases 0.8 (purple-pink) and the rest of the phases 0.3-0.7. Notice the spectral differences on the 
left-hand side during orbital minima and how the spectra are different for the rest of the phases. 



L is the luminosity of the ionizing source, n e the 
electron density of the plasma and r the distance 
between the absorber and the ionizing source. The 
fitted spectra of the orbital minima and maxima 
are presented in Figure 6, while the fitted param- 
eters are given in Table 2. 

Table 2: Spectral parameters derived 
from the fits to the orbital max- 
ima and minima. The model used is 

WARMABS* WABS*PCFABS* (MEKAL+MEKAL+GAUSS) . 

The gaussian emission line energy is fixed at 6.4 
keV. Nn of the warmabs model is the equivalent 
hydrogen column density of ionized warm absorp- 
tion. All the errors are given in 99% confidence 
level. 
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Fig. 5. — The variation of ./Vh parameter of the 
partial covering absorber PCFABS2 on top and 
the variation of the temperature parameter of the 
hardest emission component mekal3 at the bot- 
tom. Both variations are over the orbital phase. 



5 



Table 1: Spectral parameters derived from the fits to the spectra for each orbital phase of 0.1 in the 0.3-10 
keV range. All the spectra are fitted with a composite model of three collisional equilibrium plasma emission 
models at different temperatures (mekal), 1 Gaussian centered at 6.4 keV, 2 partial covering absorber model 
for intrinsic absorption (pcfabs) and a simple absorber for the interstellar absorption (wabs). The given 
errors correspond to 90% confidence level for a single parameter. 
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5. Discussion 

The orbital modulation in CVs provides es- 
sential information on the structure of the ac- 
cretion disk in non-magnetic systems and inter- 
mediate polars. The orbital-phase resolved anal- 
ysis requires demanding conditions in the X-ray 
regime such as adequate count rates, sufficient ob- 
servation times and suitable orbital/spin periods. 
Previous works on orbital phase resolved analysis 
(e.g., Pckon k Balman 2011, Hellier et al. 1996, 
Itoh ct al. 2006, Salinas & Sclegel 2004, Staude et 
al, 2008) showed that the spectral characteristics 
such as absorption column densities, plasma tem- 



perature distributions, line widths from IPs can 
be studied in detail over the orbital phase which 
contributes to the understanding of the detailed 
structure of the accretion phenomenon of these 
systems. In this work, we tried to incorporate the 
orbital-phase resolved spectroscopy in order to cal- 
culate the absorption explicitly from the disk and 
the accretion curtains/columns, separately. The 
short P sp in of the system is an advantage since 
the spectral contributions from the spin modula- 
tion can average out in a single orbital phase bin. 

It was pointed out by Evans et al. (2004) that 
the orbital modulation of FO Aqr may be a re- 



6 




0.5 1 2 5 10 ^ 0.5 1 2 5 10 

Energy (keV) Energy (keV) 



Fig. 6. — The composite spectral fit including a warm absorber fitted to the spectrum at orbital minima of 
the source (phases between 0.9 and 1.2) is shown on the left. On the right, a plot of the composite model 
alone without the data is presented. 



suit of the absorption from vertical structures on 
the disk extending out to 140° on the plane of 
the accretion disk and 25° above the accretion 
plane (however not explicitly calculated). More- 
over; Parker, Norton & Mukai (2005) showed that 
the modulation depth increases at lower energies 
which supports the absorption effects from the 
disk. The folded light curve over the orbital pe- 
riod shows three features: a plateau feature with 
almost constant /slighly increasing flux between 
phases 0.3-0.6; a peak at phases 0.6-0.8, and a 
quick decline and deep feature at the rest of the 
phases where flux minimizes around phase 0-0.1. 
In our orbital phase resolved analysis, absorption 
column of the second partial absorber (PCFABS2) 
follows this trend inversely, where it peaks during 
the orbital dipping phases and is lower during the 
orbital peak as shown in Figure 5. The absorp- 
tion column density of the maximum and mini- 
mum phases are Nn value of 17.45 xlO 22 cm -2 at 
the maximum reaching up to 5 times the minimum 
value of 3.75 xlO 22 cm" 2 . The covering fraction 
of this absorber shows no significant modulation 
over the orbital phase and is close to 1. Hence 
the absorber is almost uniformly obscuring the X- 
rays in the line of sight. The absorbing material 
is well spread out on the disk since absorbing col- 
umn never vanishes over the orbital period, but is 
denser and colder at dipping phases causing more 
absorption (See Figure 7). 

Osborne & Mukai (1989) suggest that the op- 
tical eclipse occurs at the inferior conjunction of 
the secondary star. Preceding the optical eclipse 




Fig. 7. — The schematic diagram showing the sys- 
tem. The absorber (represented as the dark region 
on the disc) is concentrated around the accretion 
impact zone (phase 0), and smeared out gradually 
around the accretion disk. 

on the binary plane by 0.2 orbital phase is a re- 
processing site noted by these authors, which is 
most likely the accretion stream-accretion disk in- 
teraction zone. When we phase-lock the X-ray 
light curve to the optical light curve using the 
cphemerides given by Osborne & Mukai (1989) 
HJD 2446081.3037 + 0.20205956E for the optical 
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eclipse, we find that the X-ray orbital minimum is 
at phase 0.8. Our phase 1.0 or 0.0 is phase-locked 
to the UV minimum given by Evans et al. (2004). 
This is the same phase 0.8 in the optical band. 

The absorption column for the simple photo- 
electric absorber component (wabs), Ah, has val- 
ues around 0.1 xlO 22 cm~ 2 and the first partial 
covering absorber component has an Ah value 
around 20xl0 22 cm~ 2 and the covering fraction 
is around 0.6-0.8 . Considering the variation of 
these three absorption components over the orbital 
phase (see Table 1), it can be inferred that each 
has different origins. The simple absorber with 
relatively lower column density than the others 
(i.e. wabs) account for the interstellar absorption. 
The first partial covering absorber (pcfabsI) with 
changing covering fraction and relatively high ab- 
sorbing column accounts for the absorption arising 
from the accretion curtain/column since it does 
not vary with the orbital phase. And finally the 
second partial covering absorber (pcfabs2) ac- 
counts for the absorption from the bulk material 
on the disk (edge) since it changes with the or- 
bital motion as explained above. At any orbital 
phase in our analysis, we can see the combined ef- 
fects of these three absorption components. The 
interstellar component is in effect at all phase in- 
tervals since it is always in the line of sight; and 
the absorption from the accretion curtain/column 
is averaged out in each orbital-phase bin since a 
phase bin of 0.1 spans an interval of 1745.8 s which 
is greater than the spin period of the white dwarf. 

The behavior of the spin pulse profile over the 
orbital period of the system also supports this sce- 
nario of the absorption components. The spin 
pulse shape have a sinusoidal nature that is con- 
stant throughout the orbital motion. The ampli- 
tude of the spin pulse vary in line with the orbital 
folded light curve. Therefore the variations over 
the spin and orbital periods have different origin, 
former by the absorption from the accretion cur- 
tains/columns and latter by absorption from the 
material on the disk near the accretion stream im- 
pact zone. 

The average spectrum of the source can be 
represented with 3 plasma emission components; 
a soft component around 0.15 keV, a medium 
component around 3 keV and a hard component 
around 30 keV. The hard component temperatures 
obtained from the orbital phase resolved spec- 



troscopy ranges between 15-70 keV. The phase 
averaged XMM-Newton data yields 13. 9±|;f keV 
(Evans et al. 2004), Suzaku data yields 14j^;° 
keV (Yuasa et al. 2010) and Swift-XRT and 
INTEGRAL-IBIS data yields 29. 7±™;* keV 
(Landi et al. 2009) which are all consistent with 
this range. The first two softer plasma temper- 
atures do not follow a significant variation over 
the orbital phase. However, during the orbital 
maximum at around phase 0.8, the flux of the 
soft component peak at about 10 times more than 
the average. On the other hand, the temperature 
of the hard component peaks during the orbital 
dipping phases. These characteristics indicate a 
signature of spectral softening during the orbital 
maximum and hardening during the orbital mini- 
mum. 

Hellicr et al. (1993) and Parker, Norton & 
Mukai (2005) argued that the orbital modulation 
in IPs may be similar to those of Low Mass X-ray 
Binaries (LMXBs) where non-symmetric struc- 
tures on the accretion disk raised by the impact 
stream (possibly around the hot spot) causing en- 
ergy dependent absorption hence reducing the X- 
ray flux. In particular a class of LMXBs called 
Low-mass X-ray binary dippers show consider- 
able fluctuations and modulation around the or- 
bital minimum. Such characteristics have been ex- 
plained by modeling the dip and non-dip spectra 
using highly ionized warm absorber models (Boirin 
et al. 2005, Diaz-Trigo et al. 2006, Balman 2009 
and the references therein). We tried to explore 
the warm absorber scenario for FO Aqr extract- 
ing spectra from phases around the orbital minima 
and maxima and fit them with the assumed com- 
posite spectral model including a warm absorber. 
We simply replaced the absorption component due 
to the disk (PCFABS2) and one of the MEKAL com- 
ponents with warmabs model. Both minima and 
the maxima could successfully be modeled with an 
ionized absorber. The ionized absoption N warma b s 
value increases by 4 times during the orbital min- 
ima which is expected since the absorber is denser 
and colder at the orbital " dip" . The values for the 
ionization parameter log(£) is consistent within er- 
ror margins between the minima and the maxima. 
This may be due to the fact that we used a range 
of 0.3 phases for the minimum and maximum cal- 
culations. The N warma bs values we find are lower 
than those found in LMXBs by about 10 times 
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at all times (min or max) (see Diaz-Trigo et al. 
2006). Since £ is given by L/(n e r 2 ), £ values 
can be compared using £ = L/NnRd where n c ~ 
Nfj/Rd (R-d, disk radius and Njj equivalent hydro- 
gen column density of ionized absorption). The 
0.6-10 keV luminosities of LMXB dippers listed 
by Diaz-Trigo et al. (2006) are of the order of 
10 36 erg s" 1 . The X-ray flux of FO Aqr in 0.6-10 
keV range is 1.4xl0~ 10 ergs cm" 2 s" 1 , yielding a 
luminosity of the order of 10 32 erg s" 1 at 400 pc 
source distance. Therefore, one can approximate 
log(£cv)^ -3 log(^LMXB) (given the warmabs ab- 
sorption column is a factor of 10 more in LMXB 
dippers). We also assumed similar P or b and Rd, 
disk radius. This yields a range of log(£cv)^ - 
0.8-1. Therefore the orbital dip in the X-ray light 
curve of FO Aqr is analogous to the dips observed 
in the dipping LMXBs with the expected cause of 
the dips being variation of temperature and den- 
sity of the warm ionized material on the disk. A 
schematic representation of the system is shown in 
Figure 7. 

6. Conclusions 

We have presented X-ray orbital-phase resolved 
analysis of the intermediate polar FO Aqr. This 
work improves upon the previous works by cal- 
culating the spectral parameters over the orbital 
phase. The distinction between the absorbing 
components are clarified and the values are explic- 
itly calculated. The absorption originating from 
the polar regions of the white dwarf can be re- 
solved from the absorption by structures on the 
accretion disk. The shape of the spin pulse profile 
is unaffected by the orbital motion, however the 
semi-amplitude of the profile change over the or- 
bital phase. The X-ray orbital variation over the 
orbit in the system arises from absorption by the 
bulge material on the disk spread well over the 
disk. The absorption column is greatly enhanced 
during the X-ray orbital dip. Moreover, we have 
modeled the absorption from the orbital dip with 
a warm absorber model for the first time for this 
source and also for CVs, confirming the ionized 
nature of the material on the disk causing the ab- 
sorption. We derived a range of ionization parame- 
ter log(£cv)^ -0.8-1 for plausible warm absorbers 
on CV disks. We calculate that the orbital modu- 
lations in FO Aqr (and plausibly some other CVs) 
are similar to those seen in LMXB dippers as vari- 



ations (density/temperature) of a warm absorber 
at the accretion impact region and spreading on 
the disk. 
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